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MODEL  FOR  ESTIMATING  METEOROLOGICAL  PROFILES 
FROM  SHIPBOARD  OBSERVATIONS 


INTRODUCTION 

Large  quantities  of  shipboard  meteorological  observations  made  throughout  the  world 
ocean  over  an  impressive  length  of  time  are  limited  in  nature  and  accuracy.  But  they  provide 
a marine  climatological  data  base  which  would  be  too  expensive  to  duplicate  to  obtain  other 
types  of  meteorological  data  subsequently  needed  but  not  specifically  included  in  the  data 
base.  One  example  of  subsequent  need  is  the  statistical  characterization  of  detailed  humidity 
and  temperature  profiles  within  the  lowest  2 kilometers  above  the  ocean  surface.  The  rela- 
tively small  group  of  radiosonde  data  taken  by  various  weather  ships  and  researchers  is  un- 
fortunately sparse  in  both  a geographical  and  a statistical  sense.  There  exists  therefore  a 
need  to  invent  a model  which  will  provide  an  estimated  profile  of  humidity  and  temperature 
when  given  the  standard  set  of  shipboard  observables.  Profiles  obtained  with  such  a model 
will  be  inferior  to  those  observed  directly  by  the  use  of  radiosondes  and  are  not  intended  to 
replace  them,  only  to  supplement  them  in  areas  where  sufficient  profile  measurements  do 
not  exist. 

The  chief  advantage  is  that  profiles  can  be  generated  where  no  profiles  have  been  mea- 
sured and  the  results  can  be  believed  to  within  a certain  degree  of  accuracy.  One  area  of 
interest  in  which  such  a model  could  be  useful  is  the  transmission  of  optical  energy  over 
various  nonhorizontal  paths.  Along  these  paths  humidity  estimates  are  valuable  both  for  its 
effect  on  molecular  absorption  of  optical  energy  and  in  its  effect  on  the  growth  of  aerosols. 

This  report  describes  such  a model  and  a method  and  data  to  test  this  model.  The  test 
method  can  also  be  used  on  future  models  to  ascertain  their  usefulness. 


THE  EMPIRICAL  MODEL 

The  model  divides  the  standard  shipboard  weather  observations  into  essentially  four 
cases: 

• fog, 

• light  precipitation, 

• heavy  precipitation,  and 

• no  precipitation. 
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Characterization  by  the  model  of  all  but  the  last  of  these  items  is  trivial  in  that  present 
values  on  the  ship  are  estimated  to  exist  at  higher  altitudes.  In  the  last  case,  however,  more 
elaborate  predictive  methods  are  used. 

An  integer  is  delivered  to  the  main  calling  program  which  gives  a rough  estimate  of  the 
quality  of  the  particular  value  being  calculated.  When  many  assumptions  are  made,  the  qual- 
ity integer  (IQ)  becomes  low  or  even  negative.  Values  of  10  are  considered  reasonable 
estimates. 

The  stability  of  the  marine  boundary  layer  is  an  important  input  to  the  model.  The 
measure  of  this  parameter  that  is  available  to  the  model  is  the  air/sea  temperature  difference. 

In  the  no-precipitation  case  the  following  distinct  classes  of  profiles  are  modeled: 

• Unstable  atmosphere  with  cloud  cover  greater  than  25%, 

• Unstable  atmosphere  with  cloud  cover  less  than  25%, 

• Neutral  atmosphere  with  cloud  cover  greater  than  25%, 

• Neutral  atmosphere  with  cloud  cover  less  than  25%, 

• Stable  atmosphere  with  a stratus  deck,  and 

• Stable  atmosphere  with  no  stratus  deck. 

The  particular  class  of  model  chosen  depends  on  the  sign  of  the  air/sea  temperature  differ- 
ence, which  determines  the  stability,  and  the  level  and/or  type  of  cloud  cover  reported.  The 
specific  characterization  of  each  of  the  profiles  from  any  of  these  classes  depends  also  on 
the  specific  values  of  the  other  shipboard  observables. 

In  the  modeling,  potential  temperature  and  mixing  ratio  are  chosen  as  the  particular 
meteorological  pair  of  variables  to  describe  the  vertical  thermal  and  vapor-loading  character- 
istics of  the  atmosphere  because  of  their  insensitivity  to  adiabatic  processes.  One  limiting 
feature  of  this  model  is  that  it  is  designed  to  describe  only  the  layer  of  the  ocean  atmosphere 
at  or  below  2 kilometers;  thus  attempts  to  use  the  model  outside  of  this  limitation  may  result 
in  severe  errors. 

The  output  of  the  model  is  relative  humidity,  absolute  humidity,  and  air  temperature 
(all  quantities  which  are  used  directly  in  propagation  calculations).  In  making  these  estimates, 
many  conversions  are  required  between  various  meteorological  variables.  Thus  a number  of 
conversion  functions  and  subroutines  are  required  by  the  model  to  convert  easily  from  one 
set  of  parameters  to  another.  These  subroutines  and  functions  (shown  in  the  listing  of  the 
program  in  the  Appendix)  use  for  the  most  part  standard  meteorological  relationships. 


Unstable  Atmosphere  With  Cloud  Cover  Less  Than  25% 

In  modeling  the  class  of  profiles  identified  as  unstable  atmosphere  with  cloud  cover  less 
than  25%  four  levels  are  considered.  Starting  from  the  layer  closest  to  the  sea,  they  are 
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• The  guperadiabatic  layer, 

• The  homogeneous  layer, 

• The  transitional  layer,  and 

• The  cloudless  cloud  layer. 


The  extent  of  each  of  these  layers  varies  according  to  the  individual  input  parameters.  The 
unstable  model  is  based  on  the  description  of  the  moist  marine  layer  in  the  trade-wind  zones 
sketched  by  Roll  [ 1 J . 

Figure  1 is  a conceptual  representation  of  profiles  of  the  potential  temperature  and 
the  water-vapor  mixing  ratio  r for  this  class  of  atmosphere.  The  first  layer  above  the  ocean 
surface  is  called  the  superadiabatic  bottom  layer.  In  this  layer  the  potential  temperature  and 
mixing  ratio  obey  a log-linear  formulation,  with  the  defining  parameters  being  related  to  the 
shipboard  observations.  The  height  of  this  layer  above  the  sea  surface  is  calculated  by  an 
empirical  formula  fitted  to  a series  of  measurements  made  by  Brocks  [2] : 


Z = 11.92  ♦ 9.69  ln(|  T - T\),  5 < T ■ ~ T , 


= T.„  - Tm.  O < Tair  - Tieg  < 5. 


where  ZJa  is  in  meters  and  Tajr  and  Tgea  are  the  air  temperature  and  sea-surface  temperature 
in  degrees  Celsius. 


The  meteorological  parameters  of  0 and  r at  any  altitude  Z within  the  superadiabatic 
layer  are  computed  from 


and 


0(Z)  = 0(0)  + 0J\nZ-^  + 4.8-fl 


Z + Zo  Z 

r(Z ) = r(0)  + rjln-g— ^ + 4.8|-] 


In  these  formulas  0(0)  and  r( 0)  are  the  values  of  potential  temperature  and  mixing  ratio 
computed  from  measurements  of  the  temperature  of  the  sea  surface.  The  variables  0,  and 
rt  are  computed  from 

0*=^W[<?(1O)  “ 0(0)1 


and 


V^io 


(r(10)  - r(0)] , 


r 


0.38 
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where 


with  0(10),  r(10),  and  u(10)  being  the  values  of  potential  temperature,  mixing  ratio,  and 
real  wind  speed  obtained  at  the  nominal  mast  height  of  10  meters.  Values  of  ZQ  are  obtained 
as  follows:  If  the 


real  wind  speed  is  measured,  then 


and  if  the  real  wind  speed  is  not  measured,  then 


Finally  the  values  of  the  so-called  stability  length  L are  calculated  from 

-r0v^“(10>2 


where  the  sea-surface  temperature  Tq  as  well  as  potential  temperatures  0^  (0)  and  0^  (10)  are 
in  absolute  temperature  units.  If  values  of  the  wind  speed  sure  not  measured,  L is  set  at  —100 
meters.  When  the  relative  humidity  is  to  be  estimated  at  an  altitude  less  than  Zn , then  it  is 
csdculated  by  use  of  0 and  r from  these  equations.  If  Z is  greater  than  Z , the  values  of 
potential  temperature  0 and  mixing  ratio  r at  Z are  used  to  estimate  tfiese  same  quantities 


The  height  of  the  lifting  condensation  level  ZQQn  must  next  be  calculated  from  measure 
ments  available  on  the  ship.  The  method  used  in  this  model  is  based  on  the  solution  of  the 
equation  for  adiabatic  lifting: 


where  T is  the  temperature  calculated  at  the  top  of  the  superadiabatic  bottom  layer,  e is  the 
vapor  pressure  calculated  at  the  same  point,  Tc  is  the  temperature  at  the  condensation  level, 
and  egc  is  the  saturation  vapor  pressure  at  the  temperature  of  the  condensation  level.  The 
potential  temperature  at  the  top  of  the  superadiabatic  layer  is  known,  because  both  the 
temperature  and  the  altitude  are  known.  Also  known  is  that  the  potential  temperature  re- 
mains constant  with  respect  to  altitude  changes  for  adiabatic  processes.  Therefore  at  the 
condensation  level  the  pressure  pc  can  be  calculated  from 


0 = Tc(1000/pc)°  . 

The  model  then  uses  the  NACA  standard  low-level  atmosphere  to  convert  from  pressure  to 
altitude. 
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Since  the  area  is  clear  (total  cloud  amount  is  less  than  25%),  then  a transition  layer  is 
considered  to  start  at  an  altitude  of  about  80%  of  the  lifting  condensation  level.  The  transi- 
tion layer  starts  at  the  t-.i  of  the  homogeneous  layer.  In  the  homogeneous  layer  both  the 
potential  temperature  and  the  mixing  ratio  at  all  altitudes  remain  constant  at  the  values 
obtained  for  the  top  of  the  superadiabatic  layer. 

According  to  measurements  by  Malkus  [3) , the  relative  humidity  in  the  subsident  non- 
cloudy  cloud  layer  is  78%.  Therefore  in  the  transition  layer  the  model  slightly  increases 
the  potential  temperature  linearly  and  adjusts  the  mixing  ratio  so  that  78%  relative  humidity 
is  achieved  at  the  altitude  calculated  as  the  lifting  condensation  level.  Throughout  the  cloud- 
less cloud  layer  the  model  sets  the  potential  temperature  to  that  calculated  at  the  lifting 
condensation  level  and  correspondingly  decreases  the  mixing  ratio  to  keep  the  relative 
humidity  constant  at  78%. 


Unstable  Atmosphere  With  Cloud  Cover  Greater  Than  25% 

Figure  2 is  a diagram  of  the  same  unstable  atmosphere  but  with  clouds.  For  altitudes 
below  the  transition  layer  the  model  is  identical  to  that  for  the  cloudless  case.  The  model 
for  the  cloud  case,  however,  keeps  both  the  potential  temperature  and  the  mixing  ratio  con- 
stant with  respect  to  altitude  throughout  the  transition  zone.  Above  the  lifting  condensation 
level  when  the  relative  humidity  is  100%,  the  potential  temperature  remains  constant  but 
the  mixing  ratio  is  reduced  just  enough  to  keep  the  relative  humidity  at  100%. 

The  cloud  layer  extends  from  the  lifting  condensation  level  to  about  2 kilometers, 
where  the  moist  marine  layer  is  usually  capped  by  an  inversion  which  divides  dry  warm  air 
on  top  from  the  cool  moist  air  of  the  marine  layer.  Above  the  2-km  level,  temperature  de- 
creases at  the  standard  lapse  rate,  and  moisture  is  allowed  to  decrease  rapidly  toward  the 
zero  level. 


Neutral  Atmosphere 

The  neutral  atmosphere,  profiled  in  Fig.  3,  is  similar  to  the  unstable  atmosphere  in 
both  the  cloudy  and  noncloudy  conditions  except  that  no  superadiabatic  layer  exists.  The 
temperature  of  the  atmosphere  is  considered  adiabatic  from  the  sea  surface  up  to  the  lifting 
condensation  level.  The  potential  temperature  measured  at  the  ship  is  used  as  the  potential 
temperature  throughout  the  lower  atmosphere  from  the  sea  surface  to  the  lifting  condensa- 
tion level  in  the  cloudy  case  and  up  to  the  bottom  of  the  transition  zone  in  the  subsident 
cloudless  case. 

The  mixing  ratio  is  interpolated  linearly  between  values  corresponding  to  relative 
humidity  of  98%  at  the  sea  surface  and  that  measured  at  the  assumed  mast  height  of  10 
meters.  Above  10  meters  and  below  the  lifting  condensation  level  the  mixing  ratio  is  set 
constant  to  that  measured  on  the  ship. 
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Stable  Atmosphere 

The  basis  for  the  model  estimates  of  the  stable  atmosphere  is  a series  of  tethered  bal- 
loon flights  made  aboard  the  USNS  Hayes  during  the  1975  fog  cruise  [4] . If  a stratus  deck 
is  observed,  the  model  constructs  the  following  atmospheric  thermal  and  water-vapor  struc- 
ture from  four  additional  shipboard-measurable  inputs:  the  sea  surface  temperature,  the 
dew-point  temperature  and  air  temperature  at  20  meters,  and  an  estimate  of  the  cloud-base 
altitude. 

The  potential-temperature  structure  is  estimated  in  the  following  way  and  as  shown  in 
Fig.  4.  The  atmosphere  is  divided  into  four  layers  separated  by  four  reference  levels:  the  sea 
surface,  20  meters  (the  height  of  shipboard  measurements),  the  cloud-base  level  (an  input  to 
the  model),  and  the  cloud  top  (assumed  to  be  500  meters  above  the  cloud  base).  In  the 
lowest  layer  a linear  interpolation  is  assumed  between  the  potential  temperatures  of  the  sea 
surface  itself  and  the  air  at  ship  height.  In  the  second  layer  a linear  interpolation  is  made 
between  the  potential  temperature  at  the  20-meter  level  to  that  at  the  cloud  base  level, 
which  is  assumed  to  be  that  of  the  ship  plus  2 °C.  Over  the  next  500  meters  an  increase  of 
3°C  of  potential  temperature  is  assumed.  Finally  above  the  cloud  layer  the  potential  tem- 
perature is  allowed  to  increase  at  a small  rate. 

Similarly  the  mixing  ratio  is  constructed  in  a piecewise  linear  fashion  in  the  same  four 
levels.  In  the  lowest  layer  a vapor  pressure  of  98%  of  the  saturation  vapor  pressure  at  the 
sea-surface  temperature  is  assumed  at  the  sea  surface.  The  mixing  ratio  is  then  linearly  inter- 
polated between  the  mixing  ratio  r t calculated  at  the  sea  surface  and  the  mixing  ratio  r2 
calculated  at  the  mast  height.  When  a stratus  cloud  is  known  to  exist,  then  saturation  is 
known  to  be  achieved  at  the  cloud  base  and  therefore  the  mixing  ratios  r3  and  r4  which  are 
necessary  for  saturation  are  assumed  both  at  the  cloud  top  and  at  the  cloud  base.  Linear 
interpolations  are  then  made  between  these  various  fixed  points.  Above  the  cloud  top  the 
mixing  ratio  is  allowed  to  drop  off  at  a 1/Z 2 rate. 

The  stable-atmosphere  no-stratus-cloud  profiles  are  identical  to  the  stable-atmosphere 
stratus-cloud  profiles  except  that  the  mixing-ratio  profile  is  assumed  constant  from  the  ship- 
mast  height  to  500  meters  above  the  calculated  lifting  condensation  level. 


IMPLEMENTATION  OF  THE  MODEL 

The  model  is  in  the  form  of  a Fortran  function:  RELHUM  (Z,  IQ,  ISHIP,  TAIR, 
VAPLD).  The  outputs  of  the  model  are  the  relative  humidity  in  percent  which  becomes  the 
functional  value  itself,  the  air  temperature  represented  by  parameter  TAIR,  and  the  water- 
vapor  load  of  the  atmosphere  represented  by  the  variable  VAPLD.  The  integer  output  vari- 
able IQ  is  a quality  index  which  is  roughly  related  to  the  trustworthiness  of  a particular 
prediction  as  described  above.  The  model  requires  only  two  inputs:  the  altitude  Z (meters)  at 
which  the  meteorological  values  are  desired  to  be  known  and  an  integer  array  ISHIP  which 
contains  12  integers  from  the  standard  shipboard  weather  reports. 

The  elements  of  array  ISHIP  are  defined  as  follows: 

ISHIP  (1)  is  an  integer  between  0 and  99  (coded  in  accordance  with  Ref.  5)  which 
describes  the  present  weather  situation; 
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ISHIP  (2)  is  a parameter  between  0 and  10  which  describes  the  lower  cloud  types 
(ILCL); 

ISHIP  (3)  is  a signed  integer  between  —999  and  +999  which  is  10  times  the  air/sea 
temperature  difference; 

ISHIP  (4)  is  a signed  integer  between  —999  and  +999  which  is  10  times  the  air 
temperature  (IAT); 

ISHIP  (5)  is  a noncoded  integer  which  is  the  estimated  height  of  the  cloud  base  of 
the  lowest  cloud  layer  in  meters  (ICLHT). 

ISHIP  (6)  is  a signed  integer  between  —999  and  +999  which  is  10  times  the  dew- 
point temperature  at  deck  level  (IDPT); 

ISHIP  (7)  is  a signed  integer  between  —999  and  +999  which  is  10  times  the  mea- 
sured sea-surface  temperature  (ISST); 

ISHIP  (8)  is  an  integer  between  0 and  99  which  is  the  wind  speed  in  knots; 

ISHIP  (9)  is  a integer  between  90  and  99  (coded  in  accordance  with  Ref.  5) 

which  represents  the  visibility; 

ISHIP  (10)  is  an  integer  which  if  equal  to  0 indicates  that  wind  speed  was  measured; 

ISHIP  (11)  is  an  indicator  of  the  horizontal  visibility  at  ship  level,  with  0 indicating 

that  visibility  was  measured  and  1 indicating  fog; 

ISHIP  (12)  is  a integer  between  0 and  9 (coded  in  accordance  with  Ref.  5)  which 
indicates  the  total  cloud  amount  (ICA). 


TEST  OF  THE  MODEL 

A test  of  the  accuracy  of  this  model  within  the  marine  boundary  layer  has  been  devised 
by  comparing  the  model-produced  estimates  with  measured  profiles  of  temperature  and 
water  vapor.  Field  data  suitable  for  this  test  was  obtained  from  the  EOMET  cruise  of  the 
USNS  Hayes  during  May  and  June  1977.  During  this  cruise  kite-balloon  soundings  were 
made,  producing  many  observations  within  the  marine  boundary  layer.  Comparisons  be- 
tween the  measured  parameters  of  air  temperature,  relative  humidity,  and  vapor  density  and 
those  same  parameters  calculated  from  the  model  which  used  the  appropriate  standard  ship- 
board measurements  taken  on  the  ship  itself  provide  a useful  test  of  any  model  of  a 
boundary-layer  profile. 

The  locations  of  the  soundings  are  shown  in  Fig.  5 superimposed  on  the  ship  track, 
which  covered  portions  of  the  North  Atlantic  and  the  Mediterranean  Sea.  Although  the 
experiment  provided  over  1700  valid  parameter  sets  obtained  with  the  NRL  boundary-layer 
sonde  [6) , the  data  are  limited  in  duration  at  any  one  location  and  therefore  do  not  form  a 
climatologically  complete  set  of  measurements.  It  is  nevertheless  an  excellent  data  base  on 
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Fig.  5 — Cruise  plot  of  the  EOMET  cruise  of  the  USNS  Hayes  during  May  and  June  1977.  The  circles 
represent  positions  during  the  track  where  kite-balloon  profile  measurements  of  air  temperature,  relative 
humidity,  and  water-vapor  density  were  made. 


which  to  test  the  model,  because  it  does  provide  a wide  variety  of  stability  classes.  The  tem- 
perature measurements  themselves  are  produced  by  dry  and  wet  thermisters  supported  on 
the  tether  line  to  the  kite  balloon  and  ventilated  by  the  relative  wind.  A solid-state  pressure 
altimeter  provides  altitude  information.  Data  are  collected  from  these  three  channels  and 

transmitted  to  the  ship,  where  they  are  plotted  in  real  time  for  the  convenience  of  operators  . ] 

and  digitally  recorded. 

The  sensitivity  of  the  system  under  real  operational  conditions  allows  both  wet-bulb 
and  dry -bulb  temperatures  to  be  measured  to  within  an  error  band  of  ±0.1 5°C.  The  absolute 
accuracy  of  the  radiosonde  is  checked  before  and  after  every  kite-balloon  flight  by  flying  the 
measurement  package  at  mast  height  and  comparing  these  telemetered  measurements  with 
the  shipboard  measurements  at  mast  height.  Thus  the  calibration  standard  for  the  airborne 
units  is  the  shipboard  Cambridge  model  100SM  air  and  dewpoint  temperature  instrument. 

The  altimeter  is  checked  before  and  after  each  flight  by  observations  of  the  pressure  at  the 
mast  height,  which  eliminates  the  time  variations  of  the  atmospheric  pressure  from  the 
altimeter.  The  operational  accuracy  of  the  altimeter  is  ±14  meters  at  the  surface  but  decreases 
with  altitude.  One  advantage  of  the  tethered  kite  balloon  is  that  the  same  instrument  package 
is  used  many  times  and  the  calibrations  of  the  devices  can  be  checked  before  and  after  each 
flight  to  make  sure  that  changes  have  not  occurred  in  the  launch  process  or  during  the  flight 
itself. 

Part  of  the  problem  in  obtaining  an  evaluation  of  the  performance  of  the  model  is  to 
include  in  the  evaluation  some  of  the  nonhomogeneous  characteristics  of  the  atmosphere 
that  might  be  present  during  a particular  test.  Sometimes  at  sea  continuously  recorded 
shipboard  observations  diverge  significantly  even  during  the  time  it  takes  to  make  a profile, 
which  is  the  time  for  the  kite  balloon  to  climb  to  its  maximum  altitude  and  to  return  to  deck 
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level.  Consequently  plots  of  all  shipboard  meteorological  observations  were  made  6 hours 
prior  to  a flight  and  6 hours  after  landing  of  the  apparatus  to  aid  in  determining  the  stability 
of  the  shipboard  observations.  The  first  part  of  the  test  was  to  prepare  from  data  for  the 
middle  of  each  flight  a list  of  shipboard  observables  such  as  would  be  entered  into  any  ship 
weather  report.  For  every  experimental  observation  of  dry -bulb  temperature,  wet-bulb  tem- 
perature, and  altitude  the  model  also  would  estimate  air  temperature,  relative  humidity,  and 
vapor  density  at  the  same  altitude,  using  as  its  input  the  list  of  shipboard  observables.  Differ- 
ences between  measured  quantities  and  the  estimated  quantities  represent  information  on 
the  accuracy  of  the  model. 

Comparisons  between  the  model  predictions  and  measured  data  have  been  made  for 
all  of  the  flights  flown  during  the  1977  EOMET  cruise.  In  some  cases  the  predicted  profiles 
of  air  temperature,  relative  humidity,  and  water-vapor  density  were  close  to  those  measured. 
Figures  6.  7,  and  8 are  examples  of  good  correlation.  They  were  obtained  from  flight  51, 
which  took  place  in  the  Mediterranean  Sea  on  June  4, 1978.  Each  point  in  the  figures  refers 
to  a sampling  of  the  data  transmitted  continuously  to  the  shipboard  receiver  and  digitally 
recorded  every  minute  throughout  the  flight.  The  scatter  in  the  kite-balloon  data  is  partly 
due  to  the  instrumental  uncertainties  in  both  altitude  and  temperature. 

Figure  9 is  a plot  of  a particularly  impressive  and  suprising  prediction  of  a dry  layer 
over  the  Mediterranean  Sea  which  was  both  measured  by  the  boundary-layer  sonde  and  pre- 
dicted by  the  model.  Here  both  the  model  estimates  and  the  measurements  show  that  a 
sharply  defined  dry  layer  exists  above  350  meters  in  altitude. 


Fig.  6 — Comparison  of  the  air-temperature  profiles  obtained  from  flight-51 
kite-balloon  measurements  (points)  and  the  model  estimates  obtained  from 
shipboard  observations  at  the  time  of  flight  51  (line) 
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RELATIVE  HUMI  01  TV  (V.) 


Fig.  7 — Comparison  of  the  relative-humidity  calculations  obtained  from 
the  dry-blub  and  wet-buld  temperature  measurements  of  flight  51  (points) 
and  the  model  estimates  of  relative-  humidity  based  only  on  shipboard 
observations  (line) 


12  13  14  15  16  17  18 

ABSOLUTE  HUMIDITY  (g/m3) 


Fig.  8 — Comparison  of  the  absolute-humidity  calculations  obtained 
from  the  dry-bulb  and  wet-bulb  temperature  measurements  of  flight  51 
(points)  and  the  model  estimates  of  the  absolute  humidity  based  only 
on  shipboard  observations  (line) 
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Pig-  9 — Comparison  of  the  model  estimates  of  a relative-humidity  pro- 
file which  predicted  a dry  upper  air  layer  from  surface  measurements 
made  at  the  time  of  flight  42  (line)  and  the  kite-balloon  measurements  of 
relative  humidity  based  on  dry-bulb  and  wet-bulb  temperatures  (points) 


i ■ 
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Not  all  predictions  are  this  reassuring  however.  Figure  10  is  a plot  of  air -temperature 
measurement  from  flight  21  near  the  Gulf  Stream.  The  air-temperature  measurements 
(although  the  most  basic  and  simplest  of  the  measurements  to  make)  seem  to  give  absurd 
data.  The  problem  results  from  the  time  required  for  the  measurement,  that  is  for  the  kite 
balloon  to  climb  to  altitude  and  then  to  be  reeled  back  to  deck  height.  Because  the  before 
and  after  calibration  measurements  of  these  data  are  accurate,  the  conclusion  is  that  the 
atmosphere  is  changing  rapidly.  Thus  during  the  50  minutes  required  by  the  kite-balloon 
system  to  climb  to  maximum  altitude  and  then  to  be  reeled  to  the  surface  the  assumption 
of  a stationary  atmosphere  does  not  hold.  Consequently  a static  model  cannot  duplicate 
these  rapid  changes  in  air-mass  characteristics.  It  may  well  be  that  the  model  prediction  is 
closer  to  a snapshot  of  the  profile  at  one  instant  of  time  while  the  profile  is  changing  from 
one  form  to  the  other  than  can  be  obtained  from  a measuring  system  which  is  slow  relative 
to  the  time  response  of  the  nonstationary  atmosphere. 

The  model  predictions  and  the  in-situ  data  taken  over  all  of  the  flights  of  the  1977 
EOMET  cruise  at  the  positions  shown  in  figure  5 were  compared  as  follows  to  allow  the  user 
some  idea  of  the  error  bars  on  the  model  estimates.  A calculation  was  made  of  the  RMS 
error  of  the  model  predictions  at  100  altitudes  equally  spaced  below  the  maximum  altitude 
obtained  on  the  flight  relative  to  the  least-square  third-order  polynomial  fits  to  the  measure- 
ment data  at  the  same  100  altitudes.  This  analysis  shows  that  the  worst-case  air-temperature 
RMS  error  by  the  model  was  3.7° C (for  flight  3)  but  that  the  average  RMS  error  for  all 
flights  was  0.7°C.  A similar  analysis  shows  that  the  model  had  an  overall  average  RMS  error 
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Fig.  10  — Example  of  an  air-temperature  profile  (points)  (obtained  in  a rapidly 
changing  air  mass  from  a sounding  system  that  takes  about  50  minutes  to  make  a 
complete  ascent  and  descent)  that  cannot  be  used  to  represent  a stationary  pro- 
file against  which  to  test  a static  model  (line).  The  data  in  this  plot  were  ob- 
tained from  flight  21. 


in  the  relative-humidity  estimates  of  about  20%.  These  errors  however  include  not  only 
inaccuracies  of  the  model  but  all  types  of  errors  which  might  exist  in  this  experiment,  such 
as  nonstationarity  of  the  air  mass. 

The  model  test  was  also  analyzed  to  see  if  the  errors  of  the  model  depended  on  altitude. 
One  would  expect  that,  since  the  model  is  based  on  measurements  at  the  surface,  the  model 
accuracy  would  decrease  with  altitude,  so  that  the  RMS  errors  between  data  points  and 
model  predictions  would  increase  with  altitude. 

Table  1 shows  the  RMS  differences  between  data  points  at  the  various  measurement 
altitudes  and  the  model  predictions  of  the  data  at  the  same  altitudes.  These  differences 
are  averaged  over  all  of  the  flights  in  100-meter  altitude  bands.  The  RMS  errors  indeed  tend 
to  increase  with  altitude. 


! 
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Table  1 — Altitude  Dependence  of  the  RMS  Errors  of  the  Model 


Altitude 

Interval 

<m) 

Number  of 
Observations 
in  Altitude 
Interval 

RMS  Error 

Air 

Temperature 

(°C) 

Relative 

Humidity 

(%) 

Absolute 

Humidity 

(g/m3) 

0-100 

480 

1.01 

9.54 

1.06 

100-200 

296 

1.21 

9.71 

0.94 

200-300 

365 

1.24 

9.86 

1.06 

300-400 

341 

1.34 

11.93 

1.40 

400-500 

155 

2.77 

20.47 

2.15 

over  500 

81 

1.65 

16.36 

1.68 
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APPENDIX 

PROGRAM  LISTING 

FUNCTION  BLLilUB  (Z,  XJ,  I SUaP,  T Aa  H,  V APLD) 


1.10EN1XF1i.ATaCN  NAME:  httriUj 

2.  TITLE:  aEIEOBOLOCIC  Al  P60Fa4.e  til  Art  AT  ES  FaOB  SHIP  OES. 

3. AUTHCN:  STUABT  U.  G AThfl  AM 
4.0BCANXZAT1UN:  NHL 

5. DATE;  JULY  1977 

6. SOUTINES  CALLED:  AbO BT , A a! , PAaT , VAFt-n  , C 1 0 , ZZ  EBO,  uZ , OSTAK, 
THZ,THSTB,FaLNT,ZCON,  SfcTXn,  au,  a EflF,  V F , T AbS , TCENT,  LIBIT 


DIOENSaON  ISHXP  (12),IP0g(14) , A Jo a a ( 1 7)  , XL LB  (o) 
1PN=aSHTP  (1) 

ILCL«ISHIF(2) 

1ASL=ISHIP (3) 

1AT»XSU1P(4) 

!CLbT=  IS  HIP  (5) 

X DP1=I  Sii  IP  (fa) 

1SS1-ISHIP  (7) 
lNS*ISHIP (8) 
x»xs=xsaxp  (9) 

1AMS*XSHXP  (10) 

XI VIS* IS HIP  (11) 

ICA*ISHXP(U) 

OAT  A (IPO  0(1)  ,1=  1,1  4) /I  0,  11,  lx, 2b, 40, 41, 42, 43, 44. 45, 

1 46,47,43,49/ 

OAT  A (TUB  IE (X)  ,1*1, 10)/14,1fa,1o,A0,Z4,50,51,S2,53,54/ 
OAT  A (IDS  IE (X)  ,1* 11, 17) /5fa, fab, fa  7, fab, S9, bO, b 1/ 

OAT  A (1CLB  (X)  ,l=1,b)/0,  1,2, 3, 4, 5/ 


CHECK  ros  XM  DICAT A OH  S OP  FOG  ABO  aF  FOUND  00  TO  100 

IF  (XVXS.  LE.94)  00  10  100 
XP(XlVXS.Ev).  1)  00  TO  100 
00  20  J*1,14 

20  XP(XPH.EwXPOU  (J))  00  TO  100 


CHECK  PCS  INDICATIONS  OP  OBIZZLt  A hu  a t 
FOUND  OC  TO  STATErtENT  NUMBE8  300 

DO  30  J*  1, 14 

30  IP(XPH.£u.  IOBIZ(J)  ) GO  TO  300 


CHECK  FOB  CLEAN  HEATHEN (IE  No  INuacaTaON  OF  PBEC.)  AND 
IF  POUNC  GO  TO  STATEMENT  NUMSa*  dOO 

\ 

DO  40  J*  1,0 

40  XP(XPn.EW>ICLN(J))  GO  TO  oOG 


XP  THE  PBOGMAB  gets  TO  THIS  PUAN t HE  MUST  A 00 ft. 

60  CALX  AbOBT (Z,iJ, XOPT,l AT,id) 

BEL  HUB*HH 
TAlN*lCE  NT (X AT) 


n c.  r.  r,  r.  n r.  r r r.  r.  r.  r>  r,  r>  r.  n r.  n r>  r>  n r.  n n r.  n r.  r.  n 


OATHMAN 


VAPID- 2.  16  5*«£La0a*VAPPiti'Ai.«i  /(iAlR»273.  1b) 
CALI  LIMITS  (UH,  VAPLD,TA1*,1W) 

BETUBB 


THIS  PAuT  UP  TuE  PBOG  BAB  IS  Pub  IlUtS  of  POu  IB  DlCA  Ti  UN  3. 

100  «ELbUB-100. 

TaIB-TCENT  AT) 

VAPID- 21 b. 5*VAPPn(TAIB)/ (iAin^B? J.  1 j) 

110  Iw»10 

IF  (Z.gT. 400.)  iw*o 
IP(Z.GT.  1000.)  Iw-2 
fitlUBB 


this  pabt  op  tue  pbogbab  is  pub  bun  prucipitatiob  times. 

bCO  IF(ILCL.Eu.b)  GO  TO  700 
IP  (Ij.cL.  Eg.  4)  GU  TO  700 
IF  (ILCL.  EU.7)  GO  TO  700 
IF  (1ASD.LT.0)  GO  TO  800 
IP  (1ASE. Eg.O)  GO  ro  yoo 


THE  SlAbLE  ATMOSPHERE 

777  CONTINUE 

IP  (Z.I.T.  10.)  GO  TO  1000 
GO  10  1100 

THIS  PAST  OF  TtaE  PBOGBAB  IS  FOB  UnuTAouE  CASES. 
DETERMINE  THE  DEPTH  OP  THE  SUPuRAiU AuATIC  LAYER. 

H00  Iw-10 

DIF-AoS ( FLO AT (IASD)/10.) 

IF IC1P.GT.0.5)  GO  TO  850 
ESA-10. • DIP 
GO  TO  B7  0 

850  ESA-11. y2»y.«>y*ALOG(DIP) 


NEXT  DETERMINE  IP  i IS  ABOVE  Ou  oELOb  THE  TOP  OP  THu 
SUPEBADIAuATIC  LAYER. 

870  IF  IE.LE. ISA)  GO  To  2000 


FOB  ESA  < 1 < EU,  THETA  Ah D w A CONSTANT  RBI  ALTITUDE 
AMD  EWUAI  10  TUE  VALUES  DETEbBIBED  a’T  a-ZSA. 

1500  T HE  A-I  bi>  (IlBS,i.NS,IO,ISSI,l  AT,  ASA) 

wSA-UE  (1  IBS,  IBS  , Iw  , I SST, ID* I, aSa, I AT) 


FIND  TUE  L1PTIBG  CO  N DE  NS  AT  lOB  UUl  FaOB  VALUES  OP  THETA 
AND  u DEItafllNED  AT  THE  Top  Op  lUE  oJRaBAUI ADATIC  LAYER. 
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C 

ZC=ZC0N (THSA,yS A.ZSA) 


DETtNBINE  THE  CLOU 0 MOUNT  Tu  IF  IulNE  IS  A 

TRANSITION  ZONE  EXISTING. 

1810  CONTINUE 

1 1 ( ICA  .L 1. 2)  GO  TO  1600 
2H*  ZC 

IN  (2. ON.  Zd)  GO  TO  1700 
18C0  E=Vf(wSA,Z) 

TAlS=TEBPlTHSA,2)-273.  1b 
ESA1*V APPB  (T  AIB) 

VAPLD»216.5*E/(TAIfi*27  3.  I) 

BELHUH*100.  *E/ES AT 

CALI  LIB ITS  (MEL8UH AFLD.TAIH, iy) 

NET UBN 


THE  IN  CLOUD  CASE 

1700  NELHUB*100. 

TAIS=TEflP(THSA,Z)-275.  15 
VAPLD=216. 5*VAPFfl(TAlH)/(TAI»*27J.  15) 
IF12.LT. 2000.)  KETUNN 
T2000=TEHP (THS A, 2000)  - 27J.15 
UlOQ0»SflIX  (Z.T2000) 

W=Ui000*4E6/2**2 
eSAI=V APPB (TAIN) 

E»Vt(U,Z) 

BELhU(l-1  00* t/ES AT 
V APLD=  zl 6.  5*E/TEHP  (THSA.z) 

NETUNN 


THIS  PANT  OP  THE  PnOGSAB  IS  FOE  A CLOUDLESS  CLOUDLAYEB. 

1600  Zd=0.8*ZC 

IF(Z.GT.ZC)  Go  TO  Is 00 
IF (Z.LT.  2d)  GO  TO  1800 
DELZ*ZC— ZH 
XZH*TEB F (IHSA.ZH) 

T2=l2H-0.00o* (Z-ZH) 

wza*ysA 

T2c*TZd— 0.006*0 ELI 

V)2C*0. 7«*SB1XN (ZC.TZC- 27j.  15) 

THSASTH (TZ#Z) 

ySA=U2H*  (yZC-yZH) • (2-2 H) /DeLe 
GO  TO  1800 


THIS  IS  THETA  AND  h CALCULATED  FUu  * SUPEN  A LAYEB 

2000  X = IbZ(IINS,lNS,Iu,ISST,IAT,Z) 

R=yZ(liNS,i»5,IU,ISST,  IDPT, l, aAI) 
TC=1EBP(X,Z)-27J.15 


GATHMAN 


VTEST-VP  (4,2) 

VAP1ST«VAPP&  (TC) 

4 EL  HUH*  1 00.  * if  TEST/ V APT  ST 
TAIR-TC 

V AP  LO* 216. 5* VIES!/ (IAltt*2  )J. Is) 
20«22ERO  (IINS,INS,  13) 

FL*tMLNT  (14S,11US,ISST,IAX,  10) 
PCN*  (2*20 )/tL 

It  (tCN.LE.O.OJ)  Go  TO  2010 
IW*lG-2 
2010  CONTINUE 

CAi.1  LIMITS  (3ELHUH  ,¥  APLD,  I AIR,  ig) 
RETURN 


: STAoit  ATHOSPUlRE,  Z>  S Hit'  Mas!  HEj.i»HT. 

iioo  iue» 

ClbAS-FLOAI  (ICLUT) 

CTOP-C LB AS*500 

It  (2.  G £.  CTOF)  GO  TO  1 1 05> 

It  (2.G£.  CLBAS)  GO  TO  1104 

T US*T.i  (TAOS  (IAT)  ,20. )♦  (Z-20.)  *2./ (Ci.BAS-20.) 
gS-SMiXR  (20.  , T C £ NT  (IDPI)) 

uo  io  :no 

1104  TUS*TM  (1  AbS(IAT)  ,20.)»2.4  (2-cLbAS)  *.00b 
US-SM1XR (20. ,TCENT (IDPT)) 

GO  10  1110 

1105  IRS*TH  (TAOS  (IAT)  ,2 0.  ) ♦ 5*. 0000  1 * (e-CTOP) 
gS*SHIXR  (20.  ,TcENT  (IDPT)  ) • Ci'0p*Cl0P/Z**2 

1110  CONTINUE 
E*VP  luS,  2) 

TA1B*TEBP(THS,2)  -27J.15 
ESA1-VAPP4  (T AIR) 

VAPLD*21b.  5*E/TEHP  (THS  ,2) 

4 ELbUH* 1 00. *E/ES AT 

It  (REi.HUB.Lt.  100.)  RETURN 

RELNUM-100. 

V AP LU*  21 b. 5* ESAT/I EBP (1HS, 2) 

RETURN 


STABLE  ATMOSPHERE, 2<  SHIP  AAoI  HE4.GH I. 

1000  XHSP*TH (TAOS  (IAT) , 20.) 

THn=TH (TABS (ISST) , 0. ) 
WSP*SUIXR(20.,TCENT(IDPT)) 
UUT*.<*B*SHlXB(1.  ,T CENT  (1SST)) 

T H S A*  X H R ♦ 2 * (IHSP— I HM) /20. 
w3A*t)Nl>2* (gSP-gNT)/20 
IU*  10 

GO  10  IbOO 


NEUTRAL  ATMOSPHERE  CASE 


S00  Iw*10 
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it  (Z.Ll.  20.)  00  ru  978 
TdSA=ld  (TAoS  (ISd  i)  ,0.) 

WSA=S«IXH  (20.  , Tc  ENT  (ID PT>  ) 

980  Lc  = 121*  lIcENT(lAT)  -TCi NT  (1DPV)  ) 
lf(2.Lt. Zc)  ou  TO  IdOO 
Oo  10  Id  10 

978  TAlfc-fLoAT (1SS  T) / 1 0. 

wSt A-0 . 9 d*SH I A (0.,lCtNT(ASSi( > 

0 H A ST=  SB  1 X (20.  / TCLNT  (1DPT)  ) 

W»  UHAST-wSEA)  *2/20.  ♦uSEA 
E«VP  (y,Z) 

LSA1=VAPPH  (1  Aid) 

UtLhUh-100.*E/E8Ar 

VAf LD*  21  o.  8*E/  irAlit*27  J.  Id) 

BtTUSN 


THIS  PANT  Of  THE  PBOOBAB  TREATS  LiOHi  «ON  FBLEZiNO  PBEC. 

OH  DB1ZZLL.  RELATIVE  HUfllDiTY  id  nSSUUED  EVERYWHERE  TO  BE 
TriAT  Of  TUB  SHIP. 

J 00  B tL BUB  *1 00 • * VAPPB  4 TO EN  T (lDP  i)  j /V APPd (TCtNl (1 AT) ) 

00  10  110 


THIS  IS  AN  APPaOXIil  ATI  ON  FOB  Idt  cLUUdLtSS  “CLOUD  LA  YEN" 
AfTEa  NALAUS  (198d) 

1900  BELbUB=7d. 

It,*  7 

TAih»TEBP(TdSA,Z)-27J.  18 
V APLD*1btt.9*VAPPB(IAlH)/(TAlB»^/d.  18) 
aETUBN 


THIS  pAaT  OF  THE  PBOOBAB  IS  fOB  THE  STdATUS  CLOUD  CASE 

700  CONUNUE 

If (2. LI. 20.)  00  TO  1000 
CLoAS* FLOAT (iCLHT) 

C IOP*c  Lb AS*  800 

If  (Z.OE.  CTOP)  OO  TO  70  d 

1F(Z.uE.  CLdAS)  OU  TO  704 

TdS*Id  (TABS  (1AT)  ,20.  ) * (Z-,.0.)  **./  (ccbAS-20.  ) 

Pc=  lH(IAbS  HAT)  ,20.)  *2 

l C*lt<iP  (PC,CLoAS)-27  J.  18 

w b*SBIXR  (20.  , TCtNl  (I DP T)  ) 

wc=SaiXB (CLbAS , TC) 

wS«wB*  (e-20.  ) • (wL-OD)/  (CcdAS-2U.) 

00  TO  1110 

704  Id S= Id  (TAdS(lAT)  ,20.)  OUo*  (--cLuAS) 

W8* SBlXB (2 ,TEBP (THS, 2) -aJj.  ID) 

00  10  1110 

708  TdSxTd  (XAbSJlAT)  ,2 0. ) ♦ 5. ♦.  Do Ow 1* ( 2-cTOP) 

w s- sm xb (clop, T ehp (Ins ,Ciop>-2  Ij. Id) *ciop*crop/u  **2 
00  10  1110 
end 


Note:  The  Malkus  citation  on  this  page  is  Ref.  3 in  the  list  at  the  end  of  the  main  text. 
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SUBBOUIINtS  AND  FUNCTIONS  FOa  Ust  Biin  aELHUN  FUNCTION. 


SUbNOUTINE  AdoBI  (E  «Iu>  A DPI#  IAT  , a a) 

X8IS  sUuBOU  TIN*.  SldPLf  TBANdBAia  Ida  SH1PDOABD  H CASU  BEHLNTS 
TO  CTHEB  ALT1IUDES  AND  DEC  BtAStj  IW  AS  T ha  PBEDICTIONS 
diccflE  LASS  ACCU BATS. 

D EM  FT=  FLOAT  (IDPT)/10. 

T ENF=FLO AT  (1A1) / 10. 

B 8=  100  .*  VAPPB  (DE  MPT)  /VAPFa  (TaNP) 

IVi=9 

IF  (2.UT.  100.)  IJ=6 
IF  (Z.vl.  500.)  Iy=J 
IF  (Z.OI.  1000.) Iu=0 
BETUBN 
END 


FUNCTION  ALT  (P) 

THIS  FUNCTION  CONFUTES  ALIITUD&  IN  NET IBS  FBOH  PBESSUBE 
IN  NldLlbABS.  THIS  FOBNUaA  IS  A FIT  TO  THE  N.A.  C.  A. 
STANDAdD  LOMEB  ATilGSPHERc.  DATA:  sdlTdSONIAN  NET. TABLES  N63 
LIST  (1968). 

IF  (F.LE.  101  J)  OO  TO  10 

AL1'=0 

BETUBN 

10  IF  (F.LE. 958.)  O O TO  20 
ALT=9.09*  (10  1J.-P) 

BETUBN 

20  ALT- 7850  *ALOU ( 1021. j 8/ p) 

BETUBN 

END 


FUNCTION  PALT(Z) 

T8IS  IS  Id*  IN  V E NS  E OF  FUNCTION  ALT  (P)  NhEBE  ALT 1TUDE, Z 
IS  IN  BaTEBS  AND  P AlT  IS  TdE  PntSoUuE  IN  HILLIDABS  AT  i. 

IF (I.L'I. SoO.)  OO  TO  10 

PALT= 102 1.  J8*EXP (-1. 27J9E-d*Si 

BETUBN 

10  PALI=101 J.-55. *2/500. 

BETUBN 

END 


FUNCTION  VAPPB (T) 

THIS  IS  AN  APrBOIi  NATION  TO  lot  Oot'F  -oBAICH  TaHPEBATUbE, 
SATUBATeD  VAPOB  PBESSUBE  OVan  li uUid  BATEN  FOBHULA. 


The  Smithsonian  tables  cited  are  Ref.  10,  and  the  List  citation  ia  Ref.  10. 
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11  MAS  DESCEIOXD  bl  BlChAnoS  l 1*71)  AMO  NIGLE*  ( 1974). 
ill  IbiS  PUNCTIOM  T IS  Ik  OtG.  e A«0  VAPPfc  IU  He. 

14-  1.0-37J./  (Xv27J.) 

HI-  13.  j1«5 
B 4-1.97o 
B 3-0.644 5 
B 4-0.1299 
HO-  1013.  25 

V APPk-JiO*EIP  (Ml*l4-H^*T4**I— aJ*i'4**j— H4*T4**4) 

BETUHH 

END 


f UNCTION  CIO  (IMS, IMS,  lv) 

THIS  PUNCTIOM  CALCULATES  IHE  dBaG  COEF ilCIENT.  IP  MIND 
SPtED  AT  10  HETEBS  IS  HEASUB*0,  OEAcONS  SUGGESTED  fOfah 
IS  USED:  BOLL  (1965),  P.  lol.  IP  MG  MIND  IS  HEASUBEU  A 
COMSTAMT  VALUE  IS  USED. 

IP  (LIUS.  He.  — 1)  Go  TO  10 
Iy-IU-1 
C 10-2E-3 
HETOBM 

10  ClQ-1.  1E-J*0.04E-3*PLOATUMS| 

BETUHH 

END 


PUNCT10M  ZZEHO (1 IMS, IM S.xwi 

THIS  PUHCTIOM  CALCULATES  THE  O In  Abie  BOUGHMESS,  ll EHO 
IH  BET BBS  AS  A PUMCTlOk  OP  MIME  SPttD  FOLLOMIHG  CHAM  MOCK 
(1955)  AHU  EXPRESSES  FeICTiOH  VCcOwITV  IH  TEfcHS  OP  THE 
DRAG  COE EPIC IE NT  AHU  <UA SUmEO  HINu  SPtED. 

If  (I1HS.  HE.-1JGO  TO  10 
ly» lu-2 
ZZEBO-  1.  E— 5 
BET IBM 

10  ZZEBO-C10(IHS,1ImS,Iu) *3. JJ  t-4 OPLOAI (IMS) **2 
BETUHH 
END 


f UNCTION  jZ(Ii.MS/IMS,IC,ISs.,  IDFI IAT) 

THIS  P UNCTION  CALCULATES  Tus.  AXE  IMG  RATIO  AT  ALT.  Z IN  AN 
AIMCSPHE BE  M Hi Ch  JBE1S  a eOG  - clMtAE  RELATIONSHIP: 

BCLI  (1965),  P.  237. 

TO-iLOAT  (1SST) /10 
yO-Q. 96* SHI IB (0, TO) 
ySH-uSlAB (IOPT,Z,IlNS,IMS,iSST,*w) 
zo- zzemo (iims.ims, iw) 

E=  U»AO)  /ZO 


The  Richardi,  Wigley,  Roll,  and  Charnok  citation*  are  Reft.  7,  8, 1,  and  9 respectively, 
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FL-FHLNT  (XUS,XIU5,  ISST.lAT.Xwi 
w2*yO«wSk*  (A  LOG  (X)  ♦N.im/FL) 
BIT  van 
BUD 


FUNCTXON  THE  (XX  US,  XUS,  Iy,  X3»T , XA1 , L) 

THaS  FUNCTION  CALCULATES  THE  POXEnTa AL  TENPERATU  BE  AT 
ALTITUDE  X XN  AN  AIHOSPHlRl  auicd  ObEYS  A LUC  - LX  NEAR 
relatxonsuxf:  roll  (1U6^  ,p.237. 

TO* a7 J. 1 S* FLOAT (XSST)/10. 

T bO *TH  (T0,0.) 
l0=2ZERO  (1XUS,XUS,  Iy) 

A*  (Z*ZG) /ZO 

FX*  FflLUX  (XUS,1XUS, XSST,XAT, lwj 
TSlAfi-THSTR (1AT, X, HUE ,IUS, XSS I, a y) 

THZ=THO» TSTAR* (ALGG(X) ♦4.o*X/F^) 

RETURN 

BUO 


FUNCTION  wSTAu  ( IDPT ,Z , 1 XUS , la S, aSS T, I y) 

TUXS  XS  TUN  "FRICTION  HIEING  RAaAUh  USED  IN  yZ  AND 
CALCULATED  USING  THE  APPROX aHaTIOnS  OF  TUE  UULk  AERO  DY NAH1C 
UtTfcOD:  BOLL  (Hob),  P.L3X  ,2  72. 

yA*SHlXR (10. ,TCENT (IDpT)  ) 
yQ=.*d*SHi  Xh  (0.  , TC  ENT  ( XSSX'j  ) 

CA=C10 (XUS.lIUS.iy) 

FK=0.3U 

ySI Afl= Syfil  (CA) * (yA-yO) /Fa 

RETURN 

END 


FUNCTXON  TUSTR (I  AT ,Z ,A XUS, XuS, XuSi , Iy) 

THIS  IS  TuE  "FRICTION  POILaTaAL  X'EflPER  ATUKE"  USED  IN  THE 
CALCULATXON  OF  TUl  USING  Tul  APpuOAIUATIOH  OF  TUE 
toULX  AtRODYNAHIC  HEThOO:  AULL(1SbDj  P.252,272. 

TA= a7J. 1 5AFLOAX (1AI)  /I  0. 

10*273.1 5 AFLOAT (XSST) / 10. 

CA=  C 10 (XUS,lXUS,ly) 

FK*C.4 

TH:»lfi=Sy  RX  (CA)  • (TU  (TA,  10. i -XU (X0,D)  j /Fk 

RETURN 

END 


FUNCTXON  FHcNT ( XUS ,1 XU S, XSS X, X AT , iy) 

C 

C l'HXS  XS  AN  APFaOXiHATlCh  FOB  Tut  fli  A IN G LENGTH  USING 

C ROlL  (IsbS) , P.  144,252. 

C 


The  Roll  citation  it  Ref,  1. 
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if  (TINS.  NE.-1)  GO  TO  10 

r blnt*-ioo 

K.-1Q-2 

NETUBN 

T0*x7 j.»FLOAT(ISST)/10. 

TA*x7J. ♦ FLOAT (I AT)  /10. 

CA-C10 (IUS,IIHS,Iw) 

X*T0*S|jHI (CA) *0.  26  46  *F  LOAX (IaS)*«2 

11*9. 8*. 3U*  (TH (TO ,0.) -Tri (iA, 10.)) 

FBLNT*— 1 .*X/X 

NETUBN 

END 


FUNCTION  ZCON (T ASA, BSA,ZSa) 

THIS  FUNCTION  CALCULATES  TUt  LIFTING  CONDENSATION  LEVEL 
NUEBE  1ASA  IS  THE  POTENTIAL  TEBFENATUBE  IN  KElVIN  AT  ALT 
ZSA  ABO  NS  A IS  TUE  SUING  aAlxO  AT  TUIS  LEVEL  IN  G/KG. 
TuIS  FOBHULATION  IS  AN  EHFalCAL  Pi 1 TO  DATA  IN  TUE 
SBIlUSONIAN  BET.  TABLES  F.32tt. 

TSA*I£BP  (TASA,ZSA) 

ESA* VP  (BSA.ZSA) 

FK*  1./.286 

X*£XP (FK*ALOG (TS A)  ) 

TC* 7 3b. 02— 150. 41  * ( ALOG 10 (X) -ALOOlO (ESA) ) 

TC-IC*7.  21*  (ALOG  10  (X) -ALOG  1 0 (ESA) ) **2 
TC*IC*27J. 15 

PZ*1000*EXP(-FK*ALOG (TASA/Xc)  ) 

ZCOk-ALT  (PZ) 

NETUBN 

END 


FUNCTION  SB1XB(Z#T) 

CALCULATES  THE  S Al'UB  AT  ION  axIING  BATIO  IN  G/KG  AT  HEIGHT  Z 
AND  TEBPEUATUHE  T (C) . 

SBIXB=0.62b*VAPPB(T) / (BAlT ( Z) -VAPPB (T) ) 

NETUBN 

END 


FUNCTION  TU(T,Z) 

CONVEBTS  TEHPEBAIUBE  ( K)  i'0  POTENTIAL  TEflPEBATUBE  (K)  AT 
ALTITUDE  Z. 

TU*T*tXP  (0. 2Uo*ALOG ( 1000./PAcl  (Z) ) ) 

TH*IH 

NETUBN 

END 


FUNCTION  TEBP(TE,Z) 


The  Smith  Ionian  table*  cited  on  thi*  page  are  Ref.  10. 
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CCHVEBTS  POTBHTlAi.  TEBPEkAXUBE  AT  ALTITUDE  Z (fl)  TO 

TEHIEkATUBE  (A) 


TEBI-TE/  lEAP(0.2Ub*ALOG(100C./BA*.X  («)  ) )) 

aiTuaa 

xao 


puhctxoh  vpcs.Z) 

COM  VESTS  BIXIHG  SATIO  (G/BUj  AT  Ai.Ti.TUOE  Z (B)  TO 
VAPCS  PSESSUBE  (BE). 

1P>PAi.T(Z)  -k/(.b22*B) 

SET  USB 
EBO 


PUHCTXOH  TABS (I) 

COM  VESTS  AS  IHPUT  TBBPEBATU»t  IHTEGEB  TO  TEBPEBATUBt  (K) 

IAbS-27 J . 1 b*PLOAT (I)  /1 0. 

BETUBB 

BHD 


FUHCTXOH  TCLHT  (X) 

COHVEBTS  AH  XHPUT  TEflPEBATUkt  XHTtGEB  TO  TEBPESATUBE  (C) 

ICtHT-PLOAT  (X)/10. 

BETUSH 

EBO 


SUUBOUTIHE  LXBXTS (Btt,V ,T,iW) 

XP ( Sri. LE. 6b. ) GO  TO  10 
IP (SB. LB. 100.)  kETUBH 
SB- 100. 

GO  10  20 
10  BB-bb. 

20  Ic-lU-4 

V«2.1bb*SH*VAPPk  (T)/(T»27j.  lb) 

BETUBH 

BHD 
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